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ABSTRACT 



The accreting X-ray pulsar Her X-1 shows two types of long-term variations, both with periods of ~35days: 1) Turn-on cycles, a 
modulation of the flux, with a ten-day long Main-On and a five-day long Short-On, separated by two 0#'-states, and 2) a systematic 
variation in the shape of the 1.24 s pulse profile. While there is general consensus that the flux modulation is due to variable shading 
of the X-ray emitting regions on the surface of the neutron star by the precessing accretion disk, the physical reason for the variation 
in the pulse profiles has remained controversial. Following the suggestion that free precession of the neutron star may be responsible 
for the variation in the pulse profiles, we developed a physical model of strong feedback interaction between the neutron star and the 
accretion disk in order to explain the seemingly identical values for the periods of the two types of variations, which were found to 
be in basic synchronization. In a deep analysis of pulse profiles observed by several diff'erent satellites over the last three decades we 
now find that the clock behind the pulse profile variations shows exactly the same erratic behavior as the turn-on clock, even on short 
time scales (a few 35 d cycles), suggesting that there may in fact be only one 35 d clock in the system. If this is true, it raises serious 
questions with respect to the idea of free precession of the neutron star, namely how the neutron star can change its precessional period 
every few years by up to 2.5% and how the feedback can be so strong, such that these changes can be transmitted to the accretion disk 
on rather short time scales. 
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1. Introduction 



accretion, accretion disks, - stars: Her X-1, - X-rays: general, - X-rays: X-ray binary pulsars. 



El 

The binary X-ray pulsar Her X-1 shows a number of peri- 
odic modulations of its X-ray flux: the 1.24 s pulse period, the 
i 1.70 d orbital period (through eclipses and the Doppler modula- 
, tion of the pulse period), a 1.62 d dip period, and a 35 d super- 
' orbital period. The last period is observed first through an on- 
oflT cycle with a lOd Main-On and a 5d Shor t-On (Fig. [B, 
separated by two 10 d Off-states (Tan anbaum et a l. 1972), and 
second throu gh a reproduced change in the shape of the 1.24 s 
pulse profile (iTrumper et al.ll986l iDeeter et al.ll998llS cott et al.' 
I2000h . With respect to these modulations. IStaubert et al. (2009) 
argued that there may be two -35 d clocks in the system which 
are generally synchronized by strong feedback. 

The 3 5 -day modulation of the X-ray flux is generally ex- 
plained by the precession of the accretion disk, which quasi- 
regularly blocks the line of sight to the X-r ay emitting regions 
near t he magnetic poles of the neutron star (iGerend & BoyntonI 
1 19761: [Schandl & Meyer 1994). The resulting clock, however, 
is irregular, showing deviations from regularity of up to ±10 d. 
This is generally demonstrated by the (0-C) (observed-minus- 
calculated) diagram (see Fig 2), plotting the diflTerence between 
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^ A matrix representing the pulse profile template is 
only available in electronic form at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsarc . u-strasbg . fr/viz-bin/qcat? J/A+A 



Fig. 1. Average 35 d flux profile of Her X-1, generated by accu- 
mulating light curves observed by RXTE/ ASM for 35 d turn-ons 
around binary phase 0.2 (Fig. 3 of Klochkov et al. 2006). The 
vertical lines indicate the binary eclipses. 



the observed turn-ons and the calculated turn-ons, under the as- 
sumption of a consta nt period (e.g . 20.5 x Pu = 34.850 d; see 
'Staubert et al.l 119831: IStill & Bovdl 120041: IStaubert et aP l2007h . 
This diagram has very interesting properties that warrant further 
attention: 1) a quasi-periodic modulation with a period of ~ 5 yr, 
corresponding to a repeated change in the mean precessional pe- 
riod by a few percent, 2) an apparent correlation with the appear- 
ance of the Anomalous Lows, which are thought to be episodes 
of low tilt ( Staubert et al. 2006) or higher twist (Leahv & Dupui^ 
[201Q|) of the accretion disk leading to a blocking of the line of 
sight to the X-ray source for some period of time (ranging from 
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a few days to a few years); 3) a p ossible long-term modulation 
with -15.5 yr (Staubert et al.l201 lb : 4) substructure with a quasi- 
periodicity of ~ 2.5 yrs in the UHURU era; 5) a correlation with 
the X-ray luminosity (iKlochkov et al .1120091) : and most remark- 
ably, 6) an intrig uing correlation with the history of the pulse- 
period evolution (IStaubert et al .1120061) . The last two properties 
require strong physical coupling between the precession of the 
outer edge of the accretion disk and the accretion torques acting 
on the neutron star when the material at the inner edge of the ac- 
cretion disk interacts with the magnetosphere of the neutron star. 
These fac ts lend important supp ort for the feedback scenario de- 
scribed in IStaubert et al.l (2009). 

With regard to the systematic variation in shape of the 
X-ray pulse p r ofiles, the situation is considerably less clear. 
iTriimper et al.l (Il986l) (based on observations of Her X-1 by 
EXOSAT) had proposed that free precession may be the reason 
for the changing pulse profiles: due to the variation in the angle 
of the line of sight with respect to the beamed emission from the 
surface of the spinning neutron star. 

Free precession may appear as a fundamental physical prop- 
erty of rigid non- spherical spinning bodies. The simplest case 
is a spheroid with some small oblateness (a "two-axial" body) 
in which the axis of angula r velocity is not aligned w ith any 
principle axis of inertia (e.g. JSommerfeld & Klein|[r897l) . It has 
been suggested as the underlying reason for the long-period vari- 
ations, both in timing and spectral propertie s, observed in sev- 
eral neutron sta rs (Jones & Andersson 2001; ICutler et al.ll2003l: 
lLink&Epsteini r2001: Haberletal. 2006). The candidate ob- 
jects are mostly radio pulsars (including the Crab and Vela pul- 
sars) the isolated X-ray pulsar RXJ0120A - 3125 (H aberl et all 
l2006h and the accreting binary X-ray pulsar Her X-1. The ex- 
istence of free precession in neutron stars and its consequences 
for our understanding of the physics of the interior of neutron 
stars i s extensively di scussed in the lite r ature ([ Anderson & Itoh| 
1975; "Shaham" "1977"; "Alpar & Ogelman 1987; Sedrakian et al 
1999: Wass erman 2003; Levin & D' An gelo 2004; Link 2007). 
Recently, [Link ( 2007) has emphasized that the question of the 
reality of free precession in neutron stars has strong implications 
for our understanding of the properties of matter at supra-nuclear 
densities. 

The idea of free precession in Her X-1 has been taken 
up by several author s (e.g. [Sh akura et al. 1998; Ketsaris et al. 
l2000[ lPostnov"2004': 'Staubert et al.h2007i: iPostnov et al. 2012). 
Recently, Postnovetal. (2012) have successfully modeled 9- 
13keV pulse profiles observed by RXTE under a particular 
set of assumptions. However, the idea has also been ques- 
tioned on various grounds (e.g., |Bi snovatvi-Kogan e t al.lll989l: 
ISoong et al.l 119871: IScott et al.l 120001: IStaubert et al. 2010b") and 
alternative models for the generation of variable pulse profiles, 
mostly involving the inner edge of the accr etion disk and/or the 
accretion column have been proposed (e.g. iPetterson et al.]|l99ll: 
IScott et al. 2000; Leahy 2004). 

Here we present the results of a model independent study of 
the variation in pulse profiles in comparison to the variation in 
the 35 d flux modulation. We find, that the two ways to count 35 d 
cycles - observing the tum-ons and observing the change in pulse 
profile shape - track each other perfectly, following one and the 
same irregular clock. First results of this study were presented by 
[staubert etaT]l2010bl It appears, that we cannot talk about two 
diff'erent clocks anymore, but we have to conclude that there is 
really only one clock in the system. This result shines new light 
onto the question of how the pulse profiles are generated and 
what the consequences are for the concept of neutron star free 
precession in Her X-1. 
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Fig. 2. Turn-on history of Her X-1: the (O - C) diagram (us- 
ing 20.5xPorb as the period to calculate the expected tum-ons). 
Upper: The Anomalous Low around MJD 51500 (AL3) is as- 
sumed to last for eighteen 35 d cycles. Lower: Th e same assum- 
ing that there were only seventeen 35 d cycles (see IStaubert et al.l 
2009). The solid line represents the linear best fit to all data, 
defining a mean period of 34.88 d (which may be used for a 
rough ephemeris of turn-ons). 



2. 35 d flux modulation and turn-on history 

The 35 d modulation shows a sharp increase of flux towards 
the Main-On. It is called turn-on, and defines 35 d phase zero; 
it generally occurs either around bina ry phases 0.2 or 0.7 
toiacconi et al. T^3l : iLevine & Jemigan 1982; see, however, 
[Leahy & Igna 2010^ Fig. [T] shows the mean flux profile for 
phase 0.2 tum-ons as observed by RXTE/ ASM (iKlochkov et alJ 
[2QQ6h . The 35 d turn-on clock is fairly irregular, allowing the 
length of an individual cycle to be either 20.0xPorb, 20.5xPorb, 
or 21.0xPorb (Stau bert et al.i 1983) (with a small fraction of cases 
showing longer or shorter cycles). Adopting P35 = 20.5xPorb (= 
34.85 d) as the mean ephemeris period, the turn-on history can 
be described by the (O - C) diagram. Fig. 2 (upper and lower) 
shows (O - C) since the discovery of Her X-1 until today as a 
function of time0 The two versions of the diagram (upper and 

^ The turn-on times listed by iLeahv & Ignal (1201 Ql) for a restricted 
period of time are largely consistent with our values. 
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Fig. 3. Representation of the template which describes Main-On pulse profiles of Her X-1 in the 9-13 keV range as function of 35 d 
phase. Left: Pulse profiles for every 0.02 in 35 d phase. Right: 3D pulse representation for the 35 d phase range 0.00 to 0.24. The 
resolution is 128 bins in pulse phase and 48 bins in 35 d phase (resolution 0.005). 



lower) diflTer in the number of 35 d cycles assumed to have oc- 
curred during the 602 day long Anomalous Low (AL3), which 
is centered at ~MJD 51500. During this gap, the observed X- 
ray flux was very low because of the blocking by the low incli- 
nation accretion disk (optical observations indicate that X-rays 
are still generated at the neutron star surface, but not seen by 
a distant obse r ver). F or the interpretation of Fig. 2 we refer to 
IStaubert et al.l (|2009), in the following called the "Two-clocks 
paper". In Fig. 2-upper the gap corresponds to 18 cycles with a 
(short) mean duration of 19.7xPorb, which is probably the cor- 
rect physical interpretation for a continued precession of the ac- 
cretion diskEl In Fig. 2-lower the gap corresponds to 17 cycles 
with a (long) mean duration of 20.8xPorb, which can be associ- 
ated with a semi-regular, long-term clock with a mean period of 
20.58 xPorb = 34.98 d. In the "Two-clocks paper" (IStaubert et al.l 
12009^) this period is associated with an underlying clock, viewed 
to be rather regular, namely (free) precession of the neutron star, 
assumed to be responsible for the observed periodic variation in 
the shape of the pulse profiles. In this model, it is assumed, that 
the precession of the neutron star is the master clock and that the 
precession of the accretion disk is evidently locked to that of the 
neutron star by existing strong physical feedback in the system, 
(nearly) synchronizing the periods of the "two clocks". 

3. Pulse profile variations 

In observations of Her X-1 by EXOSAT in 1984. lTrtimper et al.l 
(119861) discovered that the 1.24 s X-ray pulse profiles vary 
in shape as a function of the phase of t he 35 d flux mod - 
ulation. Obse rvations by Gin^^ in 1989 (IScott et al.l l2000l) . 
by HE AP- 1 (ISoong et al.l 1 1 990l) and by RXTE starting in 
1996 (IStaubert et al.l l2010allbl) . confirmed these findings and 
added a wealth of detailed information on the combined pulse 
shape and spec tral ev olution of the pulsar's beamed emission. 
PlYtlmper et all (Il986l) had suggested that the systematic vari- 
ation in pulse shape is due to free precession of the neutron 
star: the viewing angle towards the X-ray emitting regions of 
the neut ron star varies with the phase of the neutron star pre- 
cession. IShakura et al.l (1 19981) applied a model of a precessing 
triaxial shape to pulse profiles of Her X-1 observed by HEAO- 

^ A detailed analysis of the RXTE/ ASM data of AL3 is in progress. 



i, and i Ketsaris et al.l (l2000 l) did so for profiles observed by 
RXTE/FCA. Using afl observations by RXTE from 1996 un- 
til 2005 we have verified that the shape of the pulse profiles is 
reproduced every ~ 35 days. A careful timing analysis was per- 
formed of all archived RXTE data on Her X-1 and pulse pro- 
files were generated by folding with the measured pulse periods. 
Fig. [3] shows two representations of the final template: left: sl 
set of pulse profiles for 35 d phases 0.02 to 0.26 (every 0.02), 
and right: a 3D-plot for 35 d phases 0.02 to 0.24 with a res- 
olution of 0.005 (the resolution in pulse phase is 128 bins for 
0.0 to 1.0). The variation in pulse shape is evident, most pro- 
nounced is the gradual disappearance of the left shoulder of the 
main peak. This systematic a nalysis had led to the development 
of the "Two-Clocks-Model" (IStaubert et al.ll2009l) and the suc- 
cessful modeling of the observed pulse profiles by a model of 
point- and ring-like emission from the polar caps of a neutron 
star with an ofl'set-di pole field under the assumption of free neu- 
tro n star precession (Postnovll2004l: iPostnov et al.ll2012i) . 

IStaubert et al.l (|2010lj) presented a model-independent inves- 
tigation of the periodic pulse profile variations based upon con- 
structing a template of those variations for the Her X-1 Main- 
On for photon energies 9-13 keV. A set of observations was se- 
lected, providing good coverage of the 35 d phase range -0.05 to 
0.15. This template contained flux normalized pulse profiles for 
every 0.01 in phase. Any 9-13 keV pulse profile observed dur- 
ing a Main-On could then be compared to this template and the 
35 d phase be determined i^y xl^i^ minimization). For the RXTE 
data, we found that this is generally possible to an accuracy of 
±0.02 in phase. First results of using this template with profiles 
observed by Ginga, R XTE and INTEGRAL were presented by 
IStaubert et al.l ('2010a, c), showing that the turn-on history and 
the history of variations in pulse profiles appeared to be strictly 
parallel, implying that the neutron star precession (if it is indeed 
the reason behind the pulse profile variations) is synchronized 
to the precession of the accretion disk. This raised the important 
question about how the neutron star would be able to signifi- 
cantly change its precessional period every few years by up to 
2.5% and how such changes are relayed to the accretion disk on 
rather short time scales. 

We have since refined this method by constructing an ex- 
tended template from RXTE/FCA observations of two 35 d cy- 
cles, namely cycles no. 313 (Dec 2001) and no. 323 (Nov 2002). 
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Fig. 4. Example for finding "pulse profile phase-zero" for the 
case of cycle 257 Main-On, as observed by RXTE/FCA in July 
1996. 

Together, the data from the observation of these two cycles pro- 
vide a very well sampled coverage of a c omplete Main-On in the 
35 d phase range 0.00 to 0.24 (Fig. lA.lK The details of the con- 
struction of this template and its exact definition are given in the 
Appendix. An earlier ver sion of th e template (with less resolu- 
tion) is described in Stau bert et al.N 2010b). A SD-representation 
of the new template (being a matrix of 128 bins in pulse phase 
by 121 bins in 35 d phase) is given in Fig.[3]-right. 

4. Results 

We compare the characteristics of the two 35 d modulations (of 
the flux and of the variation in pulse profiles) by plotting both 
into an (O - C) diagram (Fig.O. The underlying data are sum- 
marized in Table [T] together with the 35 d cycle number and the 
satellite which performed the corresponding Main-On observa- 
tions of Her X-1. 

For the flux modulation, we use the turn-on, also called the 
"accretion disk phase-zero" and plot the (O - C) as in Fig. 2. 
Correspondingly, we determine a "pulse profile phase-zero" for 
all the Main-Ons for which there is at least one observed pulse 
profile (for most Main-Ons there are several profiles). The "pulse 
profile phase-zero" is defined and calibrated using those cycles 
which are used to construct the template: for cycles 313 and 323 
both types of "phase-zero" are identical (by definition). For any 
other cycle "pulse profile phase-zero" is found by determining 
the 35 d phases of all available pulse profiles by comparison with 
the template and a subsequent linear extrapolation to phase zero. 
Fig.|4]shows an example for the linear extrapolation to find phase 
zero. For these extrapolations, we use values of the 35 d period 
which are correct for the time of observation, based on fitting 
se ctions of the (O - C) diagram by linear functions, as described 
in lStaubert et al.l (|201 lb : to a fairly good approximation, sections 
(of variable length) can be modeled by a respective constant pe- 
riod. Similar attem pts of straight line fittings can be found in 
[Stm&Boyd 2004 IJuruaetal.ll201l[ The typical change in the 
precession period from section to section is on the order of 1 .5%. 
At the re-appearance of Her X-1 after the Anomalous Low in 
1999/2000 (AL3), the precession continued with a rather short 
period until around MJD 52000 (see Fig. O, at which time a 
rapid increase by 2.5% occurred, on a time scale comparable 
to two 35 d cycles. We note that within the entire (O - C) data 



base, the largest change (of -4%) is found in the early data from 
UHURU, comparing the initial upward trend in (O - C) to the 
following downward trend (see Fig. 2). 

Using pulse profiles of ten other Main-On cycles of Her X-1 
observed by RXTE over the last three decades, plus one from 
Ginga, one from Beppo/SAX, two from Suzaku and one from 
INTEGRAL, we find that these systematic variations are very 
stable and reproducible. The constructed template therefore al- 
lows one to determine the 35 d phase for any observed Main-On 
pulse profile in the 9-13 keV range. 

The observational result in the form of an (O - C) diagram is 
summarized in Fig. [5] The small green points (connected by the 
solid green line) represent the observed turn-on times (equal to 
"accretion disk phase-zero"). The magenta points are the times 
of "pulse profile phase-zero" as determined from the compar- 
ison of observed pulse profiles with the pulse profile template 
and subsequent extrapolation to phase zero. For these cycles the 
corresponding green points (tum-ons) are enlarged and shown 
with uncertainties. The observational evidence is clear: within 
statistical uncertainties, the values for phase-zero as determined 
by the two diflTerent methods are identical (see the last column in 
Table [T] which gives the diflTerences between the two phase-zero 
values in units of standard deviations). This means: the ''pulse 
profile clock" is just as irregular as the ''turn-on clock". Both 
clocks appear to be perfectly synchronized. 

We emphasize here, that the above result is completely 
model independent, it is obtained using only observational data. 
However, the latest results from the continued eflfort in modeling 
the observed pulse profiles by a model assumin g free precession 
of the neutron star (iPostnovl [20041: iPostnov eTal.ll2012i) leads to 
the same conclusion. This does suggest the need to abandon the 
concept of "two clocks" and to assume the existence of just one 
underlying clock which controls the variations of both 35 d zero- 
phases: that of the turn-ons, and that of pulse profile phase zero. 

Fig. \5\ demonstrates that the Turn-On clock is fairly noisy, 
with additional quasi-periodic variations on a ~5 year time scale. 
(O - C) correlates with the appearance of the Anomalous Lows 
(AL), a nd it also s trongly correlates with the neutron star's spin 
period (IStaubert et al..2006,) . 



5. Discussion 

What does the above result mean for the concept of (free) pre- 
cession of the neutron star in Her X-1? We distinguish between 
two assumptions. First, assuming the precession of the neutron 
star does indeed exist and is responsible for the variation in pulse 
profiles and( !) for those of the turn-on times, we would then have 
to find a physical explanation for two phenomena. 

1. How can the neutron star change its precessional period ev- 
ery few years by a few percent? 

2. How does the mechanism of synchronization between the 
neutron star precession and the accretion disk precession 
work? Can the feedback in the system really be strong 
enough to slave the disk, such that it follows on very short 
time scales (a few 35 d cycles)? 

There seem to be no external forces that are strong enough 
to change the precessional period, e.g. by applying a torque to 
the principle axis of inertia. The only possible origin could be 
inside the neutron star, that is, if glitches occur or if the complex 
physics of the interior of a highly magnetized neutron star with 
its crust and liquid core did allow for time variable phenomena 



4 



Staubert et. al.: Variable pulse profiles in Her X-1 



Table 1. Details of RXTE observations of Her X-1 used for the pulse profile analysis. 



Cycle 


Satellite 


Observed 




PP phase 




Calculated 


(O-C) 


(0 - C) of 


Difi". btw. Turn-On 


no.^ 




Turn- On 


+/- 


9 

zero 


+/- 


Turn-On^ 


of turn-on 


PP phase zero 


and PP phase zero 






(MJD) 


(days) 


(MJD) 


(days) 


(MJD) 


(days) 


(days) 


(sigma) 


181 


Ginga 


47642.20 


0.20 


47642.13 


0.11 


47637.87 


4.33 


4.26 


0.31 


252 


RXTE 


50111.80 


0.40 


50111.43 


0.40 


50112.46 


-0.66 


-1.03 


0.65 


257 


RXTE 


50285.90 


0.20 


50285.46 


0.20 


50286.73 


-0.83 


-1.27 


1.56 


259 


RXTE 


50356.00 


0.60 


50354.53 


0.20 


50356.43 


-0.43 


-1.91 


2.33 


269 


RXTE 


50704.90 


0.20 


50704.27 


0.20 


50704.97 


-0.07 


-0.70 


2.24 


271 


RXTE 


50773.80 


0.30 


50771.17 


0.80 


50774.68 


-0.88 


-3.51 


3.08 


301 


Beppo/SAX 


51826.40 


0.50 


51825.83 


0.20 


51820.28 


6.12 


5.55 


1.06 


303 


RXTE 


51895.20 


0.50 


51894.58 


0.22 


51889.99 


5.21 


4.59 


1.14 


304 


RXTE 


51929.10 


0.50 


51929.43 


0.50 


51924.84 


4.26 


4.59 


-0.46 


307 


RXTE 


52032.00 


0.50 


52031.37 


0.21 


52029.40 


2.60 


1.97 


1.16 


308 


RXTE 


52067.80 


0.50 


52066.46 


0.25 


52064.25 


3.55 


2.21 


2.40 


313 


RXTE 


52243.05 


0.50 


52243.04 


0.06 


52238.52 


4.53 


4.52 


0.02 




o vnm? 
KAlrL 


jz4j4.4U 


U.jU 


CO 


U.jU 


<1A An (^A 

JZ44 / .d4 


0. /o 


D.ZJ 


U. 1 L 


320 


RXTE 


52489.20 


0.50 


52488.73 


0.50 


52482.49 


6.71 


6.23 


0.67 


323 


RXTE 


52594.66 


0.50 


52594.65 


0.02 


52587.05 


7.61 


7.60 


0.02 


324 


RXTE 


52630.50 


0.50 


52629.48 


0.50 


52621.91 


8.59 


7.57 


1.44 


340 


RXTE 


53192.50 


1.00 


53192.03 


0.80 


53179.56 


12.94 


12.47 


0.37 


343 


RXTE 


53296.80 


0.50 


53297.17 


0.80 


53284.12 


12.68 


13.04 


-0.39 


353 


Suzaku 05 


53645.30 


0.20 


53644.43 


0.80 


53632.66 


12.64 


11.77 


1.06 


358 


Suzaku 06 


53817.75 


0.20 


53816.37 


0.80 


53806.92 


10.83 


9.45 


1.67 


373 


INTEGRAL 


54341.00 


0.50 


54339.06 


1.20 


54329.73 


11.27 


9.33 


1.50 



^ The numbering of 35 d cycles follows the convention introduced bv lStaubert et al.l (11983'): (O - C) = for cycle no. 31 with turn-on 
near JD 2442410. See, however, Staubert et al. 2009 for the ambiguity of counting after the long Anomaiom Low in 1999/2000 (AL3). 
^ "PP phase zero" stands for: "pulse profile phase zero" (see text). 
^ again: 20.5xPorb is assumed as the period to calculate the expected tum-ons. 
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Fig. 5. (O - C) values in units of days for observed turn-ons (green) and for so far generated "pulse profile phase-zero" values from 
pulse profile fitting (magenta), using profiles 9-13 keV profiles for several Main-Ons. The majority of the "pulse profile phase-zero" 
values are from observations by RXTE. The values based on observations by Ginga, Beppo/SAX, Suzaku and INTEGRAL are 
marked individually. 



5 



Staubert et. al.: Variable pulse profiles in Her X-1 



(see also iTriimper et al.lll986l for discussion of possible mecha- 
nisms). For the relevance of free precession to our understand- 
ing of matter at supra-nuclear densities, see e.g., Link (2007). 
The period of free precession Fpr of a two-axial symmetric body 
with moments of inertia Ii= I2 and I3 (the precession takes place 
around I3) is related to th e spin period P by the Euler equa- 
tion (Sommerfeld & Klein 1897): P^, = P I3 / (AI cos0), with 
being the angle between I3 and the total angular momentum. 
The period of free preces sion can chan ge when either AI or 
(or both) change. .Postnov et al.l (1201 2h argue that AI is proba- 
bly constant, but that a change in by about one degree, caused 
e.g. by a sudden relaxation of stresses in the neutron star crust, 
could lead to a change in precessional period by a few percent 
(assuming is of the order 50 degree for Her X-1). 

Equally challenging is the explanation of how the neutron 
star, if it is indeed the one master 35 d clock in the system, would 
be able to transmit its precessional motion to that of the accre- 
tion disk (assuming that our long-term concept is correct, that 
the turn-ons are due to the precession of the accretion disk). We 
do indeed see strong feedback in the binary system, which has 
led Staubert et al. ( 2009) to assume that the precession of the 
neutron star may be the master clock in the system, under the 
assumption, however, that the neutron star clock would be rather 
stable and that the accretion disk would have "a life of its own" 
providing the freedom to deviate from the strict clocking of the 
master - as observed in the turn-on history. We now see, that 
the turn-ons and the pulse profiles vary in strict synchroniza- 
tion, r equiring extremely strong feedback. Again, iPostnov et al.l 
(l2012h show that - under the geometrical conditions believed to 
be realized in Her X-1 - a sufficiently strong feedback could ex- 
ist through the dynamical action of the accretion stream on the 
accretion disk, such that synchronization is enforced. Certainly, 
precession of a neutron star in an accreting binary would not 
really be "free" and quantitative model calculations of the inter- 
action between the accretion disk and the magnetized neutron 
star are clearly needed. 

Second, assuming that there is no neutron star precession, we 
would associate the precession of the accretion disk with the one 
35 d clock observed in this system. The generally accepted inter- 
pretation of the 35 d flux modulation is that an inclined / warped 
/ counter-precessing accretion disk blocks our view to the X-ray 
source for most of the cycle, leaving two stretches of ~10d and 
~5 d, the Main-On and Short-On, respectively, with a more or 
less (Short-On) clear view to the X-ray source. The disk struc- 
ture may be viewed as a continuous series of rings of increasing 
radius and inclination (from inside out), shifted against one an- 
other in azimuth (angle of lines of nodes), constituting the twist 
or warp. For details of such models a nd corresponding p arame- 
ters see e.g Schandl & Mever (1994l) : IWiiers & Pringlelfl999 ): 
IScott et al.l (l2000h : lLeahvl(l2002h . It is then natural to assume that 



the disk is also responsible for the changing pulse profile shape, 
with the inner edge, precessing with the same period as the outer 
edge ( and all other rings), being the prime candidate. Scot t et all 
(l2000h . on the basis of Ginga and early RXTE observations, had 
proposed such a model in which the changes in the shape and 
in the spectral appearance of the pulses were qualitatively ex- 
plained by a combination of occultation of the X-ray emitting 
regions by the precessing inner edge of the accretion disk and 
by changes in the accretion geometry by a varying relative ori- 
entation between the disk and the neutron star magneto sphere. 
The model by Scott et al. ( 2000) requires a very small inner ra- 
dius of the accretion disk of 20 to 40 neutron star radii, which is 
incompatible with most estimates of the magneto spheric r adius 
(see e.g. the discussion and references in lScott et al.ll2000h . The 



magnetospheric radius is usually calculated assuming a spheri- 
cally symmetric dipolar magnetic field with a surface strengths 
of a few times 10^^ Gauss (as consistently found from both 
the observed cyclotron line energy of 40keV and the interpre- 
tation of the observed spin- up/spin-down behavio r by standard 
accretion torque theory, e.g. iGhosh & Lamb|[T979l IWang|[l987l 
iPema et al ]|2006>). In addition, we like to point out that the mag- 
netospheric radius in Her X-1 should not be much smaller than 
the co-rotation radius, since the source is known to live close 
to its equilibrium with respect to spin-up/spin-down. The co- 
rotation radius is near 200 neutron star radii. 

For completnes, we mention an earlier model which tries 
to explain the variation in the pulse profile shape between 
the M ain-On and the Short-On, as proposed by iPetterson et al.l 
("1991): here a portion of the inner edge of the warped, counter- 
precessing accretion disk, rotating synchronously to the neutron 
star, is raised out of the plane of the disk due to near vertical 
magnetic pressure. In this model the change in the pulse pro- 
files results entirely from the disk motion. This is consistent with 
our observations of the variations throughout the Main-On and 
supports the idea of having only one 35 d clock, namely accre- 
tion disk precession. But again, the model is only qualitative. A 
detailed quantitative model for the generation of variable multi- 
component pulse profiles by accretion disk/magnetospheric ef- 
fects is not available so far. 

Apart from the details of how the pulse profiles are gener- 
ated, the following basic questions need to be answered: 

1. What causes the initial tilt and warp of the disk and what 
keeps it alive over long time scales? In addition to initial 
conditions (e.g. the relative orientation of neutron star spin 
and orbital angular momentum at birth of the X-ray binary), 
the existence o f internal forces on the accretion disk, like 
coronal winds (Schandl & Mever "1994), radiation pressure 
(Pringle 1996; Wijers & Pringle 1999), or the impact of the 
accretion stream (Shak ura et al.ll 19991) have been proposed to 
produce and maintain the observed configuration. 

2. What causes the variation of the period of precession (by 
a few percent) on time scales from a few years to tens of 
years? - as apparent from the (O - C) diagram (Fig. 2). The 
accretion disk is subject to a large number of torques. In 
addition to the three already mentioned above, (1) coronal 
wind, 2) radiation pressure, 3) impact of the stream, there 
are three more: 4) torque due to the interaction of the neu- 
tron star magnetosphere with the inner edge of the accretion 
disk, 5) tidal forces by the binary components, and 5) the in- 
ternal viscosity of the disk. A careful analysis of the overall 
balance of these torques is needed. But it is evident, consid- 
ering the various feedback processes working in the binary, 
that both, the whole system and the accretion disk, live in 
a "delicate equilibrium", which allow s for deviations from 
strict regularity (see the discussion in IStaubert et al.ll2Q09h . 
It may even be possible that the observed ~5 yr modulation 
(and its overtones of ~10yr and -15 yr) represent a natural 
ringing freq uency of a system o f several coupled physical 
components (IStaubert et al.ll2009b . 

3. Finally, we can ask whether a system described above could 
be secularly stable. Given that the total time span of the ex- 
isting data base is short (~ 40yrs), the observations can not 
answer the question. For time scales < 40 yrs the ((9 - C) 
diagram (in the form of Fig. 2-lower) suggests that there is 
a mean precession period which could be stable, and there 
seems to be a "back-driving force" which brings the (O - C) 
always back to the center line, possibly representing a long- 
term period. 
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There is, however, evidence that the accretion disk preces- 
sion may not be stable, or at least - should a long-term sta- 
bility indeed exist - there can be temporary deviations from 
stability. The upper part of Fig. 2 shows a dramatic step in 
(O - C) around MJD 51500 (AL3). We believe (see footnote 
2 and IStaubert et al.l l2009l) that this is the correct descrip- 
tion of the physical behavior of the accretion disk, namely 
that the precessional period of the accretion disk during AL3 
and its immediate vicinity was low. So, over the time of 
this Anomalous Low the accretion disk has performed "one 
extra precessional cycle" as compared to the more regular 
clock (F ig. 2- l ower) , which in the "Two-clocks model" of 
fStauber t et al.l (l2009l) is associated with the precession of the 
neutron star. AL3 is a very special event in the long-term 
evolution of the system and it is not clear whether the notion 
that the system found back to the usual (albeit non-perfect) 
regularity after the AL is an adequate description. A com- 
prehensive study is needed about the question of stability, 
which is beyond the scope of this contribution. If it turns out 
that the system can not be stable on long time scales and that 
an external clock is needed, we would see two possibilities: 
1) Either: free precession of the neutron star, as discussed 
earlier. This would also answer the question of how the disk 
remains in phase with its 35 d clock during long Anomalous 
Lows, where the disk lies in or close to the orbital plane. 
But one would still have to ask, how the free precession can 
be stable on long time scales, when it does vary on a time 
scale of a few years. 2) Or: an intrinsic activity cycle of the 
optical companion HZ Her, leading to variable mass trans- 
fer rates with (quasi) periodic properties (Vrti lek et al.fcOOll: 
IStill & Jurua.2006: Staubert et al.. 2006. 20091) . 

6. Conclusion 

Based on our analysis of the variations in pulse profiles of Her X- 
1 with 35 d phase, as observed for a large number of Main-Ons 
by several X-ray satellites (with RXTE being the primary source 
of data), we conclude that these variations follow the same irreg- 
ular clock as the flux turn-ons. There seems to be only one -35 d 
clock in the system which is responsible for both, the flux modu- 
lation and the variation in pulse profile shape. The two phenom- 
ena are perfectly synchronized and follow each other on fairly 
short time scales. 

The important physical question is whether the underlying 
35 d clock is free precession of the neutron star which could also 
be responsible for the pulse profile variations. In this case, it 
needs to be explained how the neutron star can change its pre- 
cessional period every few years by a few percent, and how the 
neutron star can slave the accretion disk so tightly, such that the 
flux modulation follows immediately. If there is no free preces- 
sion of the neutron star, and the clock is due to the accretion disk 
alone, then the two most important questions are how the varia- 
tions in pulse profile are generated and whether the precession of 
the accretion disk can be secularly stable without the existence 
of an outside force. 
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Fig. A.l. Main-On light curves of 35 d cycles no. 313 (Dec 01, 
red) and cycle no. 323 (Nov 02, blue) from RXTE/FCA, 9- 
13keV. The count rate in cycle 313 is ~30cts/s lower than in 
cycle 323, here the points are shifted downward by additional 
20 cts/s for separation against the data points of cycle 323. 



Fig. A.2. Selected pulse profiles of Her X-1 of 35 day cycles of 
Dec 01 (cycle no. 313) and Nov 02 (cycle no. 323) as a function 
of 35 day phase (RXTE, 9-13 keV). 



Appendix A: Generation of the pulse profile 
template 

The X-ray pulse profiles of the 1 .24 s pulsation of Her X-1 show 
systematic variations with phase of the 35 day modulation. Here 
we describe the process of generation of the pulse profile tem- 
plate for the Main-On of Her X- 1 in the 9-13 keV range, its char- 
acteristics and the way in which the template is used to find the 
35 d phase of any 9-13 keV pulse profile. An ear lier version of 
thetemplate (with less resolution) is described in lStaubert et al.l 
(l2010bl) . To construct the template we use data of the Main-Ons 
of two 35 d cycles, taken by RXTE/FCA: cycles 313 and 323 
of Dec 01 and Nov 02, respectively (for a definition of pulse 
profile cycle counting see Staubert et al., 2009). Together the 
data provide a near continuous coverage of a complete Main- 
On of Her X-1. Fig. lA. II shows the two light curves in the en- 
ergy range 3-20 keV range (adjusted to a common phase zero). 
Fig. IA.2I shows examples of pulse profiles for progressive 35 d 
phases. 

For generating the template we proceeded through the fol- 
lowing steps of analysis: 

1 . Generate pulse profiles with 128 phase bins by folding events 
in the 9-13 keV range from integration intervals of typically 
0. 1-0.2 days (~0.003-0.006 in 35 d phase) with the predeter- 
mined best pulse period. From the available data we find 28 
pulse profiles at diflTerent 35 d phases throughout the Main- 
On. 

2. Normalize all profiles such that the minimum is zero and the 
maximum is 100 cts/s. 

3. If necessary, shift profiles such that the "sharp edges" of all 
profiles coincide at pulse phase 1.0. The sharp edges are de- 
fined by the fast decrease in flux after the right shoulder of 
the main peak i nto th e eclipse-like trough (see Staubert et al., 
2009b and Fig.|A31-upper-left). 

4. Record the 35 d phase (center of observation interval) for 
each profile. As 35 d phase zero the observed tum-ons are 
used: MJD 52243.05 for cycle no. 313, and MJD 52594.66 
for cycle no. 323. 

5 . Insert all normalized/shifted profiles (each profile is one line) 
into a matrix according to the corresponding 35 d phase 
(which runs along the columns). 



6. Perform cubic best fits (see below) with the 35 d phase as 
variable (along the columns) for each of the 128 pulse profile 
bins. 

7. Calculate template values using the best fit parameters with 
a resolution of 0.002 in 35 d phase. 

Cubic fits of the observed normalized count rate values are 
performed for each of the 128 columns using the function 

F (cts/s) = A -h B X (0-C) -h D X (0-C)2 -h E x (0-C)^ 

with (p being the 35 d phase (zero phase is turn-on). Fig. IA.3I 
shows the evolution of the four fit parameters with pulse phase. 
The 35 d reference phase was fixed to C = 0.1. For comparison, 
also the pulse profile of 35 d phase 0.122 is shown. 

The final 9-13keV template is a smoothed analytical rep- 
resentation of the observed pulse profile matrix. The template 
covers the Main-On in the 35 d phase range 0.00 to 0.24, it is 
a 128 X 121 matrix of normalized count rate values: each of the 
121 lines contains a normalized pulse profile (with 128 bins) cor- 
responding to 35 d phases ranging from 0.00 to 0.24 in steps of 
0.002. A graphical 3D-representation of the 9-13 keV template 
constructed from the data of 35 d cycles 313 and 323 is shown 
in Fig.[3]-right. The template matrix is electronically provided as 
an ASCII file. 

Any 9-13 keV profile observed during a 35 d Main-On can 
then be compared to this template and the 35 d phase can be 
determined (^y xiiin fitting). For the typical RXTE pulse profile, 
we find that this is generally possible to an accuracy of ±0.02 
in phase. For the comparison, it is, of course, necessary that the 
profile which is to be tested is adjusted in phase such that the 
"sharp edge" is at phase 1.0. 

In Fig. m we give an example of how "pulse profile phase- 
zero" is found using a set of profiles of the Main-On of cycle 257 
(July 1996) as observed by RXTE/FCA: six diflTerent integration 
intervals were defined for which pulse profiles were generated. 
The normalized profiles were then compared with the analytical 
template and for each of them the 35 d phase was determined. 
These phase values are plotted against the observing time and 
a linear fit with a constant slope corresponding to the current 
35 d period value is performed. For cycle 257 the time of "pulse 
profile phase-zero" was found to be MJD 50285.46. This method 
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Fig. A.3. Best fit parameters for cubic fits to the evolution of the count rate of normalized pulse profiles for the used 128 pulse phase 
bins (red points). The reference 35 d phase is C = 0.1. Upper left: A(constant), upper right: B(linear), lower left: D(quadratic), lower 
right: E(cubic). For comparison, also the pulse profile for 35 d phase 0.122 is shown. For the parameter A, which is basically the 
profile itself, the "sharp edge" at pulse phase 1.0 is drawn. 



works with any number of observed pulse profiles per Main-On, 
even with only one profile. 

Future work will extend the template to cover also the Short- 
On of the 35 d modulation of Her X-1 and to produce equivalent 
templates for other energy ranges. 
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